Why Kohn-Sham and Hartree-Fock orbitals are very close to each other. 
Shell structure of exchange potential in atoms. 
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It is found that, in closed-Z-subshell atoms, the exact local exchange potential of the density functional theory 
is very well represented, within the region of every atomic shell, by each of the suitably shifted potentials 
obtained with the non-local Fock exchange operator acting on the Hartree-Fock (HF) orbitals that belong to 
this shell. This explains the outstanding proximity of Kohn-Sham and HF orbitals and the high quality of 
the Krieger-Li-Iafrate and localized HF (or, equivalently, common-energy-denominator) approximations to the 
exchange potential. 
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Representing the quantum state of a many-electron sys- 
tem in terms of one-electron orbitals is simple and theoret- 
ically attractive approach. Such description is realized in 
the Hartree-Fock (HF) method [1], as well as in the Kohn- 
Sham (KS) scheme of the density-functional theory (DFT) 
J2] which is an efficient and robust tool, now routinely ap- 
plied in the calculations of electronic properties of molecules 
and condensed-matters structures, even very large and com- 
plex. Though the KS scheme is formally accurate, the one- 
body KS potential contains the exchange-correlation (xc) po- 
tential v xc , whose dependence on the electron density re- 
mains unknown. It is usually treated within the local density 
or generalized-gradient approximations (LDA, GGA), despite 
the well-known shortcomings of the LDA and GGA xc po- 
tentials (especially the self-interaction errors). Some of these 
deficiencies are removed when the exact form (in terms of 
the KS occupied orbitals) is used for the exchange part E x of 
the xc energy. The exact exchange potential v x is then found 
from E x by means of the integral equation resulting from the 
optimized-effective-potential(OEP) approach J3|(a)lj,|5| or by 
using the recently developed method based on the differential 
equations for the orbital shifts la LZD • The exact potential v x 
is free from self-interaction and it has correct asymptotic de- 
pendence (— 1/r for finite systems) at large distances r from 
the system; thus, unlike the HF, LDA or GGA potentials, it 
produces correct unoccupied states. In the DFT, the approxi- 
mation, in which the exchange is included exactly but the cor- 
relation energy and potential are neglected, is known as the 
exchange-only KS scheme — it is applied in the present in- 
vestigation. The full potential v xc can also be found by means 
of the OEP approach when the DFT total energy includes, be- 
sides the exact E x , the correlation energy E c depending on 
all (occupied and unoccupied) KS orbitals. This makes such 
computation tedious, to a level undesirable in the DFT, since it 
involves calculating E c with the quantum-chemistry methods, 
like the M0ller-Plesset many-body perturbation approach. 

Defined to yield the true electron density, the KS one- 
electron orbitals have no other direct physical meaning since 
they formally refer to a fictitious system of non-interacting 
electrons. However, it is a common practice to use these or- 
bitals in calculations of various electronic properties; in doing 



so the A^-electron ground-state wave function of the phys- 
ical (interacting) system is approximated with the single de- 
terminant built of the KS orbitals. This approximate approach 
is justified by (usually) sufficient accuracy of the calculated 
quantities, which is close to, or often better than, that of the 
HF results. It seems that the success of the DFT calculations 
would not be possible if the KS determinant, though being 
formally non-physical, was not close to the HF determinant 
which, outside the DFT, is routinely used to approximate the 
wave function of the real system. Therefore, understand- 
ing this proximity is certainly very important for the funda- 
mentals of the DFT. 

Previous calculations HI @] have shown that, not only the 
whole KS and HF determinants, but also the individual oc- 
cupied KS and HF orbitals o <t(i") and 4>™(r) in atoms are 
so close to each other that they are virtually indistinguishable 
(here the orbitals, dependent on the electron position r and 
the spin a ={, f, are numbered with index a = 1, . . . , N a ; 
Ni + 2Vf = N). This proximity is particularly remarkable 
for the exchange-only KS orbitals. It is surprising in view of 
the obvious difference between the exchange operators in the 
KS and HF one-electron hamiltonians (see below) and the fact 
that the corresponding KS and HF atomic orbital energies, e aa 
and e™, differ substantially, up to several hartrees for core or- 
bitals in atoms like Ar, Cu J4j,|5[]. This apparent contradiction 
has not yet been resolved; in Ref. 11 1 OTI it is suggested that 
the KS and HF determinants are close to each other "since the 
kinetic energy is much greater than the magnitude of the ex- 
change energy". The present paper investigates the proximity 
of the KS and HF orbitals and reveals that it results from the 
specific properties of the HF orbital exchange potentials y[] 
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obtained with the non-local Fock exchange operator built 
of the HF orbitals. Simultaneously, it becomes clear why, 
in atoms, the DFT KS exact exchange potential v xa (r) has 
the characteristic structure of a piecewise-like function where 
each part spans over the region of an atomic shell and it 
has distinctively different slope dv xa (r)/dr in consecutive 
shells [11]. These shell constituents of v xa (r) are found in 
the present Letter to be very well represented, within the re- 
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gion of each atomic shell, by the suitably shifted potentials the term (defined using the sign convention of Refs. |0,0]) 



r (r), Eq. <Q}, for the HF orbitals (j>™{r) that belong to 
this shell. Thus, it is shown how these orbital-specific HF 
potentials, which describe the exchange in the physical sys- 
tem, are almost directly mapped onto the KS local exchange 
potential v xa (r) which is common to all electrons with spin 
a in the KS non-interacting system. The revealed properties 
of v^. a (y) are also found to stand behind the high quality of 
the Krieger-Li-Iafrate(KLI) [ 3] approximation to the exact po- 
tential v xa (r), as well as, of the localized HF (LHF~) lll2ll ap- 
proximation, equivalent to the common-energy-denominator 
approximation (CEDA) II 1311 . Note that the KLI potential can 
be derived by assuming that the exchange-only KS and HF or- 
bitals are identical while the LHF (CEDA) potential is found 
when the KS and HF determinants are assumed to be equal. 

The HF one-electron spin-orbitals <fi™(r) are obtained by 
minimizing the mean value (\I r |ii"|\f') where H is the Hamil- 
tonian of the TV-electron interacting system and ^ belongs to 
the subspace of normalized TV-electron wave functions that 
are single determinants built of one-electron orbitals. Sim- 
ilar minimization is carried out in the (exchange-only) OEP 
method, but there are two additional constraints (cr =t, I) 
that for every trial determinant all N a constituent spin-orbitals 
4>aa (r) satisfy the KS equation with some local KS potential 



r (r). The minimizing potential V SCT (r) 
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(r), yields, 



after subtracting from it the external u e xt(r) and electrostatic 



s (r) terms, the exact exchange potential v xa (r) 
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(r). 



Thus, the KS equation, satisfied by the corresponding (OEP) 
orbitals </w(r) and their energies e aa , takes the form 



+ Wes "tot + V x 



(2) 



The total electron density n tot (r) = nt(r) + ni(r), which 
enters w es [n tot ] (r) = J dr' n to t(r')/l r ' — r l' i s tne sum °f tne 
spin-projected densities n a (r) = Yl a =i l</w( r )| 2 - In the HF 
equations, satisfied by the orbitals (p™(r) and energies e™, 
the multiplicative local exchange potential v xa (r) is replaced 
with the non-local Fock exchange integral operator {^(r), 
built of {<?^} =i, while the potential v es (r) is found for the 
HF total electron density defined in a similar way as rt tot (r). 
Real KS and HF orbitals are used throughout this paper, de- 
pendence on a is suppressed hereafter. Both the KS and HF 
equations need to be solved selfconsistently. 

The exact exchange potential v x = v® EP satisfies the OEP 
equation 



N 

E 

a=l 



<f> a (r)S<p a (r)=0, Vr 



(3) 



which results from the OEP minimization and depends on v x 
through the orbital shifts (OS) 5(/) a (r). Each OS fulfills the 
equation |g|7|] (below (j> a , e a are the solutions of Eq. (O) 



(4) 



and it is subject to the constraint (<fi a \S4>a) = 0. The equation 
(|4|i includes the KS Hamiltonian /i s , present in Eq. (|2), and 



Wt{v) = [vl{v) + D aa - v x (r)] a (r) 



(5) 



Here, we have v ¥ x = i>l [{&}], D aa = 



The OS 8(j) a and the constant D aa give, within the pertur- 
bation theory, the first-order approximations to the differences 
-{4>a F ~ M and -(e™ - e a ). Here, the orbitals <j>™ and the 
corresponding energies e™, are the solutions of the HF-like 
equation which is the same as Eq. (ffl except for v x replaced 
by v ¥ built of the KS orbitals {4>b}- The corresponding per- 
turbation is then equal to 6h s = v x — v x so that the first-order 
correction to <j) a {r), i.e., -5c/> a (r) = - X^i^a c «a <M r )' 
c ta = (<l>t\Sh s \<i> a ) / (e t - e a ), satisfies Eq. © and the con- 
straint indeed. Obviously, the solutions <j>™, e™ are not iden- 
tical to the self-consistent HF orbitals <f>™ and orbital energies 
e™. However, the relations ||A0 a - {-5<f> a )\\ < 0.13|| A^J 
(where ||0|| 2 = / dr |^(r)| 2 ), |Ae a - (-5e a )\ < 0.003|Ae a | 
11411 found for Be and Ar, validate the representations of 



K by 



b a and Ae Q = 



e a by 



8e a = —D aa used in the presented argument; the above in- 
equalities are obtained for <j) a , e a , 84> a calculated as in Ref. 
I2D, and 4>™ (expanded in the Slater-type-orbital basis), 
taken from Ref. fiUl . 

The part of W a {v) = Sh s (r)(j> a (r) parallel to the orbital 
(j) a is W a (r) = —D aa (j) a (r) and it sets the energy shift Se a . 
The part W^{v) = W a {r) - Wi\r), perpendicular to <p a , 
sets the OS £0 a (r), Eq. ©. Thus, the KS and HF orbitals, 
(/> Q (r), ^ F (r), can be close to each other, even if the orbital 
energies e a , e™, differ significantly, provided the term (r) 
is sufficiently small. Note that the orbitals remain unchanged 
when a (possibly orbital-dependent) constant is added to the 
Hamiltonian in the KS or HF equations. 

For closed-Z-subshell atom, the non-local (integral) Fock 
exchange operator, acting on an atomic orbital a (r) = 
r~ 1 Xni{r)Y lm (9,<f) (a = nlm), yields u£(r)0 a (r) = 
r~ 1 F x . nl (r) Y lm (9, 4>) where Y lm (8,Lp) is the spherical har- 
monic; hereafter, n, I, m denote the principal, orbital, 
and magnetic quantum numbers. The factor F x - n i(r) is 
defined |l|] with the exchange integrals and determines 
the v xa (r) = v x - n i(r) = F x -„i(r) / 'xm{r) . The OS 
64> a (r) = r~ 1 Sxni(r)Yi m (6,(p) depends on W^(r) = 
r ~ 1 ^' md ( r )Yi m (0, if) through its radial part 



F x - n i(r) + [D nl . <n i - v x (r)] Xni(r) (6) 



entering the equation for Sxni (r) found from Eq. 

Neither the OEP minimization nor the resulting OEP equa- 
tion ([3]l imply readily that the orbital shifts 5(f> a (r) are small. 
However, the occupied KS and HF orbitals differ so minutely 
from each other that, for atoms, the OEP total energy is only 
several mhartrees higher than the HF energy 75hf 
The proximity of individual HF and KS orbitals can be quan- 
tified with the norms ||<50 Q || which are indeed very small, in 
comparison with \\<fi a \\ = \\4>^ F \\ = 1- Calculating the OS 5(f> a 
with the method of Ref. [J, we obtain \\8d>u\\ = 0.00669, 
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FIG. 1: The OS norm square 



(grey bars) and the contribu- 



tions c nH . nl (stacked bars) to it from bound states 



for the 



occupied states cj>„i m in the Ar atom; the Is bars are magnified by 
the factor 20. 



0.00630 for Be and ||<5</> Q || < 0.007 for each oc- 
cupied orbital in the Ar atom. The partition ||<50 n z m || 2 = 
Y^'^n c n'ini> plotted for Ar in Fig. Q] shows that, among the 
KS bound orbitals 4> n 'hn, the dominating contributions (? n , lrd 
to the nlrn OS come from the n'lm orbitals, n' — n— 1, n+1, 
i.e., from the neighboring electronic shells; e.g., for Scf>s s in 
Ar, the largest terms c^, ; . ni are found for the n'l — 2s (occu- 
pied) and n'l — 4s (unoccupied) orbitals. But, there remains 
a large part of \\d(j) n i m || 2 which cannot be attributed to higher 
unoccupied bound states 4> n 'im since the corresponding c 2 /;.^; 
terms vanish rapidly; see Fig. Q] This unaccounted part comes 
from continuum KS states (e„/; > 0). These results confirm 
that the assumption 6<fi a = and 6(f> a — Y^t^a c ta4>t = 0, 
used to derive the KLI and LHF (CEDA) approximation, re- 
spectively, are very accurate but not exact. 
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FIG. 2: (a,b) The potentials v° EF (solid line), v x:n i (dashed and dot- 
ted lines) in the Ar atom; the HF potentials 5™ ( follow w x; „; (within 
the figure resolution); (c) The differences v° EV — sf" (dotted line), 



v °i* - „™ (dashed line), v x 0EP 



v™ OEP = « x os (solid line); the 
up and down arrows mark the points r™ n+1 and r™ m , respectively. 



The norms 



the shells S„ 

_L:rad/ 



> n, where both the functions Xni(r), 



have such low values because the W ,' m (r) decay exponentially. 



terms W^' rad (r) are sufficiently small for all r. This, com- 
bined with the relation 

v X ;nl(r) + D nknl = v x (r) + aL-U , (7) 

Xni{r) 

found with Eq. ((6), implies that each shifted KS-OEP orbital 
exchange potential 



Since the KS orbitals ^ = (j) a [ v x \ found for the 
exact exchange potential v x EP are very close to (j>™, the 
terms F x;n/ (r), v x - n i(r), and D nlM [v x ] (for any v x ) obtained 
with {<f>a EP } are virtually indistinguishable from the quanti- 
ties F ^ni( r )> w "n(( r )' D rt,nli v x\ calculated with the HF or- 
bitals {0™}. Thus, the terms W n ,' rad [u x , {0a F }](?") are very 
close to W^l' tad [v x , {</>a EP }]( r )- As a result, they are small 



ii{ r ) = v x -ni(r) + D„i ;n i 



(8) 



for v x = 



,,OEP 



and, also, by continuity, for any potential 



is very close to the exact exchange potential v x (r) = v x EF (r) 
within the r-interval (r„_i in , r n n +i) = S n where the de- 
nominators in the r.h.s. of Eq. (0, i.e., the orbitals x n ;( r ) 
from the ?i-th atomic shell (K, L, M, . . .), have largest mag- 
nitudes; the shell border points r n , n +i for n = 0, 1, . . . , n occ 
(the respective HF points r™ n+1 , defined precisely below, can 
be used) are near the positions 7-™ ln where the radial electron 
density p(r) has local minima. The potentials v x - n i(r) are also 
close to v x EF (r) in the central parts of shells S n i, n' < n, 
though not so tightly as for the shell S n . The proximity of 
the potentials, evident in Fig. |2fa,b) for the Ar atom, holds as 
well for other closed-Z-subshell atoms. It is disturbed in the 
vicinity of the nodes of Xni(r), where the potential v x , n i(r) 
diverges while the term W n ,' mi (r) is finite and small. The po- 
tential v x - n i(r) also differs significantly from v® EP (r) within 



v x close to v x E¥ . Such a class Vo [17] of potentials v x that 
yield small W^' rad [i> x , {0™}] and, consequently, lead to the 
KS orbitals </> Q [i> x ] almost identical to cf>™ [18], exists pro- 
vided it is possible, for given v x , to find constants C n i that 
make terms U n i(r) = F™ t + C n /x™ - v x x™ small for all 
r. Indeed, we then obtain D™. nl w C n i so that the term 

W ni™ d [ v *' {<^a F }] ~ U nl is small. Low magnitude of U n i(r) 
implies that, within each occupied shell S n , the shifted HF po- 
tentials v™ nl {r) = vlf nl {r) + C nl (I e C n = {0, . . . , Zi" ax }) 
lie very close to v x (r), and, as a result, they almost coincide 
with each other, 



£iW~&W, l,l'GC n ,reS n 



(9) 



(similarly, as the OEP potentials v x;n i(r) do; Fig. |2fa,b)). A 
generalization of Eq. (O is found when, in the expression for 
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U n i(r), the potential v x (r) is replaced by v™ n , v (r) for r G £„/ 
and the smallness of /7 n ; (r) is accounted for. The generalized 
relation reads 

FZl + CnlX™ « (C + X™ , re S n , (10) 

and it is satisfied for suitable set of constants {C n i} for all in- 
dices (nl), {n'V) corresponding to the occupied HF orbitals. 
This is an intrinsic property of the HF orbitals (and the Fock 
operator), since it is not implied by the DFT or the definition 
of v x EP , though it is revealed here by inspecting the results 
for v x = u° EP . The total energy E[v x ] = (^[v x }\H\^[v x }), 
for any v x € Vo, is very close to Euf (due to 4> a [v x ] ~ tp™), 
and, as a result, also to i?[t>° EP ] since the potential v x EP min- 
imizes the functional E[v x ] > Euf- This can explain why the 
exact exchange potential v x EP belongs to the class Vo and, as 
a result, it gives the KS orbitals <f> a very close to (f>™. 

Assuming that the constants C n i satisfying the relation (TlOb 
are known, we can construct a continuous piecewise potential 

occ 

<{r)=Y,0™{r)v ( ?\r) (11) 

n 

formed from the HF shell exchange potentials 

v { x n) {r) = ffi(r) + Cru] P™i{r)/pf{r) , (12) 

each applied in its shell region S n . The points r™ n+1 defining 
the shell borders are the solutions of the continuity equation 

vi n \r) = v x n+1 \r) forn = l,2,...,n occ - 1; = 0. 
We denote p«f = 4tt(2/ + l)( X ^f) 2 , P™ = ZieC n P™> 
P HF = ET Pn- = 6{r - r^jeKl^ ' r). Each 

shell potential v x n ^ (r) is very close to the almost coinciding 
potentials v^ nl {r) = v™ nl (r) + C nU I G £„, for r G S n . 

Thus, the relation ([Tol l with uSL// replaced by v x n ' means 
that the terms Uroir) are small for v x = v x w . This also holds 
for the KLl-like [3] potential 

occ 

"KLI-HF/ \ _ ( F HF _L r< ^/ HF \ ^ HF / (\ X\ 

v x ( r ) = l^«;ni + C "'X„; )X„i/ P , (13) 

ni 

since, due to Eq. ([Tol l, it is very close to v^,,, (r) w v x n ^ (r) 
for r G S 1 ,!', and, thus, to v x w for all r. Therefore, both po- 
tentials v x w and v^ Ll - liF belong to Vo- In particular, it can 
be shown to hold for v^ 1 ™ = wf J " HF [{C^ LI - HF }] where 
the constants are found from their self-consistency condition 
= A^b*™ {X?JA} given in Ref. &. 
The potential «x W (r), built of the HF shell potentials vi 71 ^ (r) 
with C„; = C^ I " HF , is found (see Fig. [2c)) to be a very 
good approximation to the exact exchange potential v x EP (r). 
The quality of its approximation is almost the same as of the 
potentials tiK LI HF and w^ LI " 0EP . The latter is defined similarly 
as vf Ll ' liF but with the HF orbitals {x™ } replaced with the 



OEP ones { X ^ ? } and with C nl = D nl , nl K° EP , {x°f P }]- The 
potential w^ LI " 0EP is the dominant part of the exact exchange 
potential v x EP — vl^ 1 ' ^ + u°f , where the minor part v x s , 
depends linearly on {5<fi a } Ja |7[]. It is this OS term v x s (r) 
that produces the characteristic bumps of v x EP (r) at the shell 
borders; cf. Fig. |2|c). 

In summary, we find that when, for each HF orbital, a suit- 
ably chosen (orbital-specific) constant is added to the Fock 
exchange operator in the HF equation, the electrons occupy- 
ing different HF orbitals are subject to very similar local ex- 
change potentials (and hence they move in very similar total 
potentials) within the atomic regions where the orbital radial 
probability densities are substantial. This proximity is par- 
ticularly tight for the exchange potentials of the orbitals that 
belong to the same shell and it holds in the region of this shell. 
Thus, each HF orbital is only very slightly disturbed when the 
(shifted) orbital exchange potential is replaced in the HF equa- 
tion with a KS exchange potential (common to all orbitals) 
that, within each shell, lies very close to this orbital potential. 
As a result, in each shell, the KS exact exchange potential is 
very well approximated with the shifted HF orbital exchange 
potentials of this shell, and, even better, with their weighted 
average - the shell exchange potential, Eq. (fT2l ). 

Discussions with A. Holas are gratefully acknowledged. 
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